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Above the saturation field, geometrically frustrated quantum antiferromagnets have dis- 
persionless low-energy branches of excitations corresponding to localized spin-flip modes. 
Transition into a partially magnetized state occurs via condensation of an infinite number 
of degrees of freedom. The ground state below the phase transition is a magnon crystal, 
which breaks only translational symmetry and preserves spin-rotations about the field di- 
rection. We give a detailed review of recent works on physics of such phase transitions and 
present further theoretical developments. Specifically, the low-energy degrees of freedom of a 
spin- 1/2 kagome antiferromagnet are mapped to a hard hexagon gas on a triangular lattice. 
Such a mapping allows to obtain a quantitative description of the magnetothermodynam- 
ics of a quantum kagome antiferromagnet from the exact solution for a hard hexagon gas. 
In particular, we find the exact critical behavior at the transition into a magnon crystal 
state, the universal value of the entropy at the saturation field, and the position of peaks in 
temperature- and field-dependence of the specific heat. Analogous mapping is presented for 
the sawtooth chain, which is mapped onto a model of classical hard dimers on a chain. The 
finite macroscopic entropies of geometrically frustrated magnets at the saturation field lead 
to a large magnetocaloric effect. 

§1. Introduction 

Theoretical investigation of finite field properties of geometrically frustrated 
magnets is a new rapidly developing direction of research, which is to a large extent 
stimulated by the experimental studies. There are now a few examples of frus- 
trated magnetic compounds such as pyrochlore Gd2Ti2C>7, 1 ) garnet Gd3GasOi2, 2 ^ 
and spinel CdCr204, 3 ) for which full or partial reconstruction of the phase diagram 
in magnetic field has been done. From theoretical point of view, an applied mag- 
netic field competes with antiparallel alignment of spins favored by antiferromagnetic 
exchange interaction and creates an extra source of frustration. This may lead to 
various interesting phenomena at high magnetic fields such as exotic quantum ground 
states, novel types of the critical behavior, the magnetization jumps and plateaus 
and so on. Here we give a detailed overview and present further development of the 
recent theoretical works on a universal behavior of geometrically frustrated magnets 
in the vicinity of the saturation field Hg.^ 10 ^ 

At high fields Heisenberg antiferromagnets exhibit a phase transition between 
a fully polarized state and a state with reduced magnetization. A partially magne- 
tized state typically breaks spin-rotational symmetry about the field direction and 
has a long-range order of transverse spin components. Since the quantum ground 
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state at H > H s is known exactly one can also calculate an exact spectrum of single 
magnon excitations u;(k). n ) For ordinary nonfrustrated antiferromagnets w(k) has 
a well denned minimum at a certain wave- vector Q with a gap A = H — H s . The 
high-field transition can be described as a Bose condensation of spin-flips with mo- 
menta k = Q. 12 ) The above simple picture fails, however, for frustrated magnets. 
Reflecting degeneracy of classical ground states at H < H s , the magnon dispersion 
at H > H s has a continuous set of degenerate minima. For two well-known frus- 
trated models, Heisenberg antiferromagnets on a face centered cubic lattice and on 
a frustrated square lattice, single magnon spectra in the saturated phase are degen- 
erate on lines in the Brillouin zone. The phase transition at H = H s is, therefore, 
associated with condensation of an infinite number of soft modes. This yields an 
effective dimensionality reduction and change in the critical behavior. 13 ) Geometri- 
cally frustrated antiferromagnets on kagome, checkerboard, and pyrochlore lattices 
exhibit even stronger degeneracy with a completely flat lowest branch of magnons: 
uj(k) = H — H s . The theory of such phase transitions is discussed below on two 
examples of the sawtooth chain and the kagome antiferromagnet, Fig. 1. 

§2. Localized magnons 

We consider nearest-neighbor Heisenberg antiferromagnets in an external field 
H = ^Ji i S i -S i -H-^S i (2-1) 

(ij) * 

with a general value S of on-site spins. The exchange coupling constants Jij connect 
only the nearest neighbor sites. The couplings are all the same for a kagome lattice 
Jij = J, while the sawtooth chain is described by two exchange constants: J\ for 
base- vertex bonds and J2 for base-base bonds, see Fig. ^ Above the saturation field 
H s the ground state of (|2-lj) is a ferromagnetic vacuum |0) = | TtT •••), where all 
spins are in a state with the maximum possible value of S 1 ? = S. The low-lying 
excitations are spin flips \i) = 5"" 1 0) . Diagonalizing the Heisenberg Hamiltonian 
(|2-1|) within one-magnon subspace, which is possible since S£ ot = Y2i Sf commutes 
with 7Y, one finds an exact one-particle spectrum. n ) 

The primitive unit cell of a kagome lattice contains three sites. Accordingly, 
there exist three branches of magnons for every wave-vector: 

ut(k)=H-6JS, W2, 3 (k) = H-3JS± JSyfifl + 2 7k ) , (2-2) 

where 7k = l/6^e lk ' ai is a sum over six nearest-neighbor sites on a triangular 
Bravais lattice. The lowest energy branch wi(k) has no dispersion and its energy 
vanishes at the saturation field H s = 6JS. Among the two dispersive branches the 
lower one uisfe) has also a vanishing gap at the saturation field. 

In the case of the sawtooth chain, one, generally, finds two dispersive branches 
of one-magnon excitations: 



u 1>2 {k) = H — 2J\S — — cos k) ± Sy Jf (1 — cos k) 2 + 2J 2 (1 + cos k) . (2-3) 




However, for a special choice of the coupling constants J2 = 0.5Ji (J\ = J), the 
lower branch becomes dispersionless: 

u 1 (k) = H-US, U2(k) = H — JS(1 — cos k) . (24) 

The above two branches are separated by a finite gap Au = 2JS. 

The two formulated spin models are related to several magnetic materials. A 
kagome lattice model is applicable, for example, for SrCrg p Gai2-9pOig, 14 ) while the 
sawtooth chain describes magnetic delafossite YCUO2.5. 15 ) Unfortunately, both 
compounds have rather large exchange constants J ~ 500-fT, which makes impossible 
to reach their high- field regimes H ~ H s . 

The dispersionless excitations correspond in real space to localized states. Spin 
flips are completely localized excitations in Ising models, while for Heisenberg interac- 
tion they hop between adjacent lattice sites and generally acquire a finite dispersion. 
In the two considered examples localization of spin flips is determined by lattice 
topology. For a kagome lattice a simplest localized state corresponds to a magnon 
trapped on a hexagon void. Its wave function is given by 

i^) = 7=fE(- 1 ) ri " 1 ^i°)' ( 2 - 5 ) 

* n=l 

where numbering of sites inside the i th hexagon goes counterclockwise starting with 
the lower left corner, see Fig. ^ An outside spin on one of the triangles surrounding 
hexagon is connected to two sites, where a localized magnon has equal amplitudes 
and opposite phases. This produces destructive interference and the amplitude of 
spin flip on surrounding sites vanishes. The energy of the localized magnon (|2-5j) is 
H — 6JS, that is the same as wi(k). More complicated eigenstates with the same 
energy are constructed by taking closed loops with even number of spins, which 
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contain only two sites of any crossed triangle. The wave-function of a localized state 
has equal amplitudes on all chosen sites with alternating phases and ir. An example 
of a spin flip localized on two hexagons is shown in Fig. ^ The wave-function of such 
a state is a linear combination of simple hexagon states: \ipj j+ x ) oc (\<fij) + \tpj+ x )). 

The localized modes on smallest hexagons can be used as a nonorthogonal basis 
in the subspace of one-magnon states from the dispersionless branch (|2-2[) . To show 
this we note, first, that there are as many hexagons Nh = N/3 on an iV-site lattice as 
the number of states in the Brillouin zone. Second, there is only one linear relation 
Yli^Wi) = between hexagon states for a cluster with periodic boundary condi- 
tions. It is constructed using the following arguments. An arbitrary site is shared 
between only two hexagons i and j. Wave- functions of the corresponding localized 
states have opposite signs for spin-flip amplitudes on the chosen site. Accordingly, 
the two amplitudes dj and a,j from a linear relation must be equal. This prop- 
erty, when extended by induction, leads to a unique linear relation between hexagon 
states: a sum of all states with equal amplitudes <Zj = const. The torus topology of 
a periodic cluster allows presence of a few other one-magnon states with the same 
energy H — 6JS, which cannot be decomposed into a hexagon basis. Such localized 
magnon states correspond to closed lines with nontrivial winding around the cluster. 
The localized magnon loops obtained as linear combinations of hexagon states have, 
in contrast, zero winding and are contractable into a point. In the thermodynamic 
limit N — > oo presence of a few states with nontrivial winding is unimportant and 
the localized hexagon states can be used as a real-space basis for the lowest branch 
in Eq. (f2~2j). 

For the sawtooth chain a localized magnon is trapped in a valley between two 
triangles. Its wave- function is 



The numbering of spins inside one valley is shown in Fig.^with identity Sn = S3 
The probability for a spin- flip to jump to an outside site vanishes in the state Ij2-fi|) due 
to the relation J\ = U<i- The N/2 valley states again form a complete nonorthogonal 
basis in the subspace of the lowest branch u\{k) (|2-4jl . Since bottom (base) spins are 
not shared between adjacent valleys, the linear independence of states Ij2-fij) is fulfilled 
even for finite periodic chains. A more extended localized state on a sawtooth chain 
shown in Fig. ^ is a linear combination of states in the two adjacent valleys. 

Localized one-magnon states allow to construct a class of exact multiparticle 
states. Configurations, where localized spin flips occupy isolated clusters, become 
exact multimagnon states of the quantum spin Hamiltonian. Note, that for nonfrus- 
trated quantum magnets in two or three dimensions, it is not possible to construct 
exact multiparticle states beyond the two- magnon subsector. 11 ) Among all exact 
eigenstates with noninteracting localized magnons, states with the highest density 
of spin flips play a special role. These are constructed as close-packed structures of 
magnons localized on smallest clusters: hexagons or valleys. For a kagome lattice the 
close-packed structure is constructed by putting localized magnons on every third 
hexagon, see Fig. |5J Such a magnon crystal with -/V max = N/9 localized spin-flips 
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Fig. 2. Magnon crystals built from localized spin flips for a kagome lattice (top) and the sawtooth 
chain (bottom). 

breaks a three-fold translational symmetry of the magnetic Hamiltonian. Transverse 
components of spin operators on sites belonging to different localized magnons re- 
main uncorrelated. The magnon crystal does not break, therefore, spin-rotational 
symmetry about z-axis. A magnon crystal for the sawtooth chain has N max = N/4 
localized spin-flips and is two-fold degenerate. 

Since isolated localized magnons do not interact with each other, it is plausible 
that they correspond to the lowest energy states in every n-magnon subsector with 
n < iV max . The proof of the above property is easily formulated for lattices with a 
uniform exchange constant Jij = J like kagome, checkerboard, and pyrochlore anti- 
ferromagnets. 5 )~ 7 ) Let us consider, for simplicity, only S = 1/2 case. The Heisenberg 
Hamiltonian (|2T|) can be split into an Ising and a transverse parts according to 

H« = jY t SfSj-HY;si i H ± = ij^(5/-57+5-5+) . (2-7) 

(u) { (ij) 

Below we always subtract energy of the ferromagnetic state from 7i zz . A quantum 
state of n isolated localized magnons is an eigenstate of 7i with E n = n(H — H s ), 
where H s = 3 J for a spin-1/2 antiferromagnet on a kagome lattice. In addition, such 
a state is also an eigenstate of H zz and H 1 - with E zz = n(H — 2 J) and = —nJ, 
respectively. The idea of the proof is to show that for an arbitrary n-magnon state: 
(H zz ) > E zz and (H^ > E±. 

For H zz the formulated relation holds trivially: once all spin flips occupy different 
bonds the Ising part of the Hamiltonian has the minimal energy: E^ z = n(H — ^zJ), 
where z = 4 is the number of nearest neighbors for a kagome lattice. In order to 
show that a similar inequality holds for the transverse part "H we map the n-magnon 
subspace of a spin-1/2 model onto the Hilbert space of n hard-core bosons B n . The 
transverse part of the Heisenberg Hamiltonian is, then, the kinetic energy of bosons 
Jt^ = ft. In addition, we define the Hilbert space of n bosons without the hard-core 
constraint B n Q with an evident relation B n C B n Q. The minimum of the kinetic 
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energy is easily found in B n Q\ all n particles must occupy one of the lowest energy 
single-particle states with the kinetic energy E^ = — J found from Eq. (|2-2|) . If the 
hard-core constraint is imposed, then, from the variational principle, the expectation 
value of K, can only increase: mia(K)B n > min(^)g n0 = — nJ = E^. 

The above proof does not immediately apply to the sawtooth chain, which has 
nonequivalent bonds. A localized magnon IJ2-6JI is an eigenstate of the total Hamil- 
tonian Id, but not of Id 1 - or 7i zz . In order to solve this problem we rewrite the Ising 
part in Eq. (|2-7|) and split the Heisenberg Hamiltonian as Ji = 7i\ + H.2, 



ij , 



- E Ji i (2 ~ S i 

(ij) 



2 j 



(2-8) 



Localized states (|2-(i|) are eigenstates of Ti.i with the eigenvalue H—2J (S = 1/2) and, 
trivially, of H.2 with zero eigenvalue. Then, the proof for the sawtooth chain copies 
the above arguments for the kagome antiferromagnet applied to 7i\ and TL2 instead 
of Ti. 1 - and Ti, zz . Similar conclusion has been earlier obtained using a somewhat 
different approach in Ref. 6). 

After establishing that isolated localized magnons are the lowest energy states, 
the high-field magnetization process of geometrically frustrated spin systems at T = 
can be understood as follows. Above H s all magnetization subsectors are separated 
from the ferromagnetic ground state by finite gaps E n = n(H — H s ). At the satu- 
ration field all gaps in subsectors with n < JV max vanish and a huge degeneracy of 
the quantum ground state develops for frustrated antiferromagnets. Below H s the 
subsector with the largest possible density of localized magnons n — A^ raax corre- 
sponds to the lowest energy. Therefore, there is a finite jump of the magnetization 
at H = H s between a saturated phase and a magnon crystal state. Such universal 
jumps with AM = N max /N have been found in numerical exact diagonalization 
studies of finite clusters of the sawtooth chain, kagome and checkerboard lattice an- 
tiferromagnets and several other frustrated one-dimensional models. 4 - 1 ' 5 )' 8 )~ 10 ) It is 
intuitively clear that noninteracting localized magnons being the lowest energy state 
in the corresponding magnetization subsectors make also the largest contribution to 
the partition function at low temperatures. Moreover, since the localized excitations 
have no dynamics their statistics is equivalent to statistics of specially chosen classi- 
cal particles. Using such a mapping one can find an exact low-temperature behavior 
of quantum frustrated magnets in the vicinity of the saturation field. 7 > Corrections 
to such an asymptotic behavior should be exponentially small once the isolated lo- 
calized magnons are separated from from higher-energy propagating states by finite 
gaps. Below we consider separately the low-temperature properties of the sawtooth 
chain and kagome antiferromagnet. 
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§3. The sawtooth chain 



Contribution of degenerate levels to the partition function is determined by di- 
mensionality of the corresponding subspace of the total Hilbert space of a quantum 
system. In order to calculate the dimensionality we can utilize an arbitrary basis 
for degenerate states. Such a basis has to be complete but not necessarily orthogo- 
nal. Among different configurations of noninteracting localized magnons only states 
formed by magnons localized in isolated valleys have to be counted for the sawtooth 
chain. Extended localized states, such as a state on the right-hand-side of Fig. ^ are 
represented as linear combinations of the basis valley states. The multimagnon 'val- 
ley' states can be associated with dimer coverings of a base chain: once a localized 
spin-flip (|2-6j) occupies a given valley or a site of the base chain it is impossible to 
put another localized magnon to the same or the adjacent valleys. The hard-dimer 
states in a given magnetization subsector are linearly independent. This property 
can be proven by induction using linear independence of the valley states (|2-6|) in 
the one-magnon subsector. In the following, we estimate gaps which separate the 
hard-dimer states from a continuum of the scattering states and then calculate the 
low-temperature thermodynamics of the sawtooth chain. 

3.1. Two-magnon bound states 

The scattering states of two magnons can be in principle calculated exactly 
using a general method discussed in Ref. 11). We shall adopt instead a variational 
approach, which is much simpler for the present problem. The two branches of single 
magnon excitations for the sawtooth chain Eq. (|2-4|) are separated by a finite gap 
Au = 2JS. The low-energy scattering states are, therefore, formed by magnons 
from the lowest branch u>i(k). In the two-magnon sector there are ^(N/2)(N/2 + 1) 
states constructed from N/2 one-magnon states of the dispersionless branch. The 
number of the basis hard-dimer states is ^(N/2)(N/2 — 3). The difference between 
the two yields the number of low-energy scattering states, which is equal to N. In 
real space representation the above scattering states correspond to localized magnons 
occupying either the same |</?f) or the adjacent \(pi(fi + i) valleys. We again assume 
5 = 1/2 for simplicity and construct accordingly the following variational basis for 
the low-energy interacting two-particle states: 



The exclusion principle of spin-flips has been used for normalization in the above 
equation. States and have nonvanishing overlaps: 
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The nonzero matrix elements of the Hamiltonian between these states are 



Z4 ~ ~ ~ 

(ipi\H\ipi) = 2H — —J , (iPi\H\^i) = 2H - 2J , (^i±i\H\^i) 



2H-AJ 
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(i>i±i\nm = \j , mn\i>i) = (^+i|w|^> = WH o * 4J . (3-3) 



A simple variational ansatz for a propagating state of two localized magnons in 
adjacent valleys is constructed as \k) = ^ e ~ lkr i It has the following energy 

v ' (k\k) y ' 35 + 2 cos fc v ; 

Two interacting localized magnons cannot scatter on each other because of an ab- 
sence of propagating one-particle states at low energies. Instead they form a bound 
two-magnon pair with energy E(k). Finite dispersion of such pairs comes from as- 
sisted hopping of adjacent localized magnons. The minimal value of E(k) is reached 
at k = with the gap A = 20/37J ~ 0.54J separating the low-energy boundary 
E2 = 2{H — 2J) from the higher energy states. Numerical minimization of E(k) 
for an improved variational ansatz \k) = e~ lkVi (\ipi) + c|V>«)) yields A ~ 0.44J 
with an optimal value c ~ —0.24. This value of the gap compares very well with 
the numerical diagonalization result 0.42J. 9 ) The gap is determined mostly by a 
nearest-neighbor repulsion of localized magnons in adjacent valleys. The dispersion 
of bound two-magnon complexes E(k) is weak and does not exceed 7%. We con- 
jecture that the intervalley repulsion leads to a similar behavior in all n-magnon 
sectors: the lowest energy states with E n = n(H — 2 J) are separated from the scat- 
tering states by a finite gap A n = O(J). At temperatures T <C min n {Z\ ra } one can 
neglect the higher energy states and consider only the contribution of isolated local- 
ized magnons. The latter problem is equivalent to a one-dimensional lattice gas of 
particles with energies H — 2 J and on-site and nearest-neighbor exclusion principle, 
which is also known as a classical hard-dimer model. 

3.2. Low-temperature behavior 

We consider classical hard dimers on a periodic chain of length L = N/2, N 
being the number of spins in the sawtooth chain. The partition function of this 
model is 

E 

T 



X^exp 

W} 



- a^) , (3-5) 

(ij) 



where \i = H s — H is the chemical potential and Oi = 0, 1 are the occupation numbers 
for dimers and fee = 1- The nearest-neighbor exclusion principle is imposed by the 
last term. The partition function 1)3-5(1 can be exactly calculated using the standard 
transfer matrix approach. In the thermodynamic limit the free energy is determined 
by the largest eigenvalue of the transfer matrix: 



^ /A r = _IrinQ + ^i + e^j. (3-6) 

The entropy S = —dJ-jdT depends on external magnetic field and temperature via 
li/T: 

S,N = i ln(l + VTTJ^) - i ln2 - £ (l - 7 ==) ■ (3-7) 
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Fig. 3. Temperature dependence of the entropy of the sawtooth chain in magnetic field. Symbols 
are exact diagonalization data, 9 ' dashed lines correspond to the hard-dimer model, solid lines 
are obtained in the effective Ising model. 

During an adiabatic process the combination {H — H s )/T remains constant and, 
consequently, temperature drops to zero as H — > H s . Such a very strong magne- 
tocaloric effect is a direct consequence of the condensation of a macroscopic number 
of soft modes at H = H S . 1%S} At H = H s the entropy remains finite down to T = 
and has a universal value: 



The localized magnons of the sawtooth chain behave, therefore, similar to param- 
agnetic degrees of freedom in an effective field h = H — H s . The major difference 
with an ideal paramagnet is that the entropy of the latter at h = depends on a 
spin length, whereas the residual entropy of frustrated magnets at H = H s has a 
geometric origin. 

Using the transfer matrix technique, one can also calculate the dimer-dimer 
correlation function. The result is 



where the average dimer density a and the correlation length £ are given by 




(3-8) 



(a i(Tj ) = a 2 + (-iy- j a(l - a)e^ n ' r ^, 



(3-9) 



a = sin 2 — = 



tan 7 = 2e^ /2T 



2 l + 4e^/ T + Vl + 4e^/ T ' 



7=2 log ( cot 




(3-10) 



The correlation length exponentially diverges £ ~ e M//2T as T — > 0. 
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Fig. 4. Temperature dependence of the specific heat of the sawtooth chain in magnetic field. 
Symbols are exact diagonalization data, 9 - 1 dashed lines correspond to the hard-dimer model, 
solid lines are obtained in the effective fsing model. 



The specific heat is obtained by taking temperature derivative of the internal 
energy £ = —1/2 [io: 

C/N = — sin 7cos7=— — T - ,„, o;o . (3-11) 

' 8T 2 ' ' 2T 2 (1 +4e^/ T ) 3 / 2 

The specific heat as a function of x = ^l/T has two peaks at x = x± determined by 
solutions of the following equation: 

x + 2 = 2e x {x-4) for x + = 4.0526 , x_ = -2.8159 . (3-12) 

Comparison of the entropy and the specific heat of the sawtooth chain calculated 
in the hard-dimer model and by a full numerical diagonalization of a cluster with 
./V = 20 sites 9 ) is shown in Figs. |2H3 The analytic and the numerical results agree 
well with each other below T* ~ 0.1 J. In particular, the hard-dimer representation 
correctly reproduces the residual entropy at H = H s (Fig. 01 . At T > 0.1J a thermal 
contribution of the higher energy states becomes significant. Those states can be 
partially taken into account by using analogy with Ising models near the saturation 
field, which are also mapped to statistical models of hard-core objects. 17 ^' For the 
sawtooth chain such a generalized Ising representation can be constructed because of 
a weak dispersion of bound two magnon pairs l|3-4jl . Once the /c-dependence of E(k) 
is ignored, there are no zero-point fluctuations and we end up again with a model of 
classical hard-core particles, where the nearest-neighbor exclusion is replaced by an 
Ising type repulsion: 

H = ^{-fJ-cfi + Vo-jCTi+i) , (Ji = 0, 1 . (3-13) 

i 

For the strength of the Ising interaction we have to choose the gap, which separates 
bound two-magnon pairs from noninteracting localized magnons: V = A ~ 0.44J. 
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The classical Ising model (|3-13p is solved exactly by the transfer-matrix method with 
the following result for the free energy: 

jr/ N = -I T InQ [1 + e- {V -^' T ] + yj^ [l - e~^-^/ T ] 2 + e^A . (3-14) 

The above expression transforms into the hard-dimer result (|3-6J) for T <C V. Com- 
parison between numerical exact diagonalization results and the two analytic ap- 
proximations is presented in Fig. |31 The entropy calculated in the Ising model shows 
a nice agreement with the numerical data up to a remarkably high T* ~ 0.3J. Note, 
that the Ising repulsion V is not a fitting parameter, but has been calculated above. 

The specific heat of the sawtooth chain as a function of temperature in a constant 
magnetic field H ~ H s (Fig.EJ) has a peak at T m = (H s —H)/x±, where ± correspond 
to H < H s and H > H s , respectively. Such a low-temperature peak appears due 
to a freezing of localized degrees of freedom and is well described by both analytic 
approaches. At higher temperatures T ~ J, the specific heat has a second peak, 
which reflects development of one-dimensional short-range spin correlations. Its 
presence is not captured, of course, by the low-temperature hard-dimer mapping. 
Experimental observation of the double-peak structure in temperature dependence 
of the specific heat of real frustrated magnetic compounds can be a clear sign of 
localized magnons. 

Field dependence of the specific heat is presented in Fig. |SJ C(H) has a deep 
minimum at the saturation field. The minimum occurs due to a weak temperature 
dependence of S(T, H s ), which is completely T-independent in the hard-dimer model. 
According to Eq. ()3-12j) . the high- and low-field peaks in the specific heat correspond 
to H m = H s —Tx±. Comparison between the three types of the results for the specific 
heat shown in Fig. |S] allows also to determine the field range, where the localized 
magnons play a dominant role: H > 0.8H S . 
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Fig. 6. Magnetization curves of the sawtooth chain obtained in the effective Ising model. Full 
lines correspond from top to bottom to T = 0.02 J, 0.05 J, 0.1J, and 0.2 J; dashed line is the 
zero-temperature magnetization jump. 



The zero-temperature jump between a half-magnetization plateau and the fully 
saturated phase is smeared at finite temperatures. According to the hard-dimer 
representation, the magnetization curves M(H,T) = (1 — <r)/2 cross at H = H s 
with the universal value M(T, H s )/N = (1 + v / 5)/4y / 5. The finite-temperature 
width of the jump is approximately AH ~ AT. The above value of the magnetization 
corresponds to an average density n = (y/E-l)/WE « 0.1382 of localized magnons. 
Magnetization curves of the sawtooth chain are illustrated in Fig. H3 by using the 
effective Ising model, which give more accurate results for the selected temperatures 
than the hard-dimer mapping. 



Following our analysis of the sawtooth chain, we define a class of hard-hexagon 
states for a kagome antiferromagnet near the saturation. They are built from 
magnons localized on smallest hexagon voids of a kagome lattice, such that no two 
hexagons have common sites. The hard-hexagon states are linearly independent 
in the thermodynamic limit with at most one linear relation in every magnetization 
subsector. The centers of hexagons form a triangular lattice and the above states are, 
therefore, mapped onto a classical gas of hard-core particles on a triangular lattice 
with the nearest-neighbor exclusion principle. The latter is a famous exactly solvable 
model 18 ) and, therefore, a great deal of the exact information is known about such 
states. However, in order to apply the available results to a kagome antiferromagnet 
we first have to show that (i) the hard-hexagon states are separated by finite gaps 
from the higher-energy states and (ii) there are no other low-energy states in every 
magnetization subsector with n < N max or their contribution is vanishingly small. 



§4. Kagome antiferromagnet 
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4.1. Two-magnon bound states 

Low energy part of the two-magnon subsector of a kagome antiferromagnet con- 
sists of ^(N/3)(N/3+l) states constructed from one-particle states of the flat branch 
u>\(k). These include ^(N/3)(N/3 — 7) hard-hexagon states. Additional 4N/3 in- 
teracting or scattering states correspond to N states of localized magnons, which 
occupy adjacent hexagons, \(pnpi +ai ), and to N/3 states with two spin flips on the 
same hexagon, \(ff). In a simple treatment of a spin-1/2 model we neglect the latter 
states and define the following real-space basis for the former ./V states: 



35 , , 

n,n'=l 

where I = 1-3 and ai = (1,0) and a2 5 3 = (±1/2, Vo/2) are three nearest-neighbor 
sites on a triangular lattice of hexagons. The normalization factor in Eq. (|4-1|) takes 
into account the exclusion principle of spin flips. The overlap matrix elements for 
the above states are 

(^2#h) = -1/7 , (^3#li) = &U±a lt \if>u) = 1/35 . (4-2) 
The Hamiltonian (|2Tj) has the following nonzero matrix elements 

(Ai\n\^u) = 2H-^-j, (foilHWii) = - 10H ~ 29J , mn\^ t ) = 2H ~ 6J . 

7 35 35 

(4-3) 

All other nonvanishing matrix elements are obtained from Eqs. (|4-2jl and 1)4-3(1 by ap- 
plying symmetry operations of a kagome lattice. Using ansatz |k) = e _ikr *q|V>zi) 
for a bound state of two magnons we calculate its energy as 

£(k) = = 2(H - 3 J) + e(k) , (4-4) 



eCk 



2 |ci| 2 + |c 2 | 2 + |c 3 | 2 -^P k 



7 (|ci| 2 + |c 2 | 2 + |c 3 | 2 )[l + ^(cosA;i + cos/c2 + cosA;3)]-ip k + ^Q k ' 

P k = Cl c* 2 (l + e~ iks ) + c 2 c* 3 (l + e ifcl ) + c 3 cl(e- ik i + e ik *) + c. c. , 

Q k = Cl c* 2 (e ik i + e~ lk *) + c 2 c* 3 (e lk * + e" ifc3 ) + c 3 c\{l + e^-**)) + c . c. , 

where k{ = k • a^. The minimum of e(k) occurs at zero total momentum k = for 
c\ = C2 = C3 with a gap Z\ = 0.24J. This value of the gap can be again used as an 
estimate for the repulsion between localized magnons on adjacent hexagons, though 
dispersion of bound pairs ~ 12% is somewhat larger in two dimensions. 

Among various states of n isolated localized magnons there are also states, which 
do not obey the hard-hexagon rules. The simplest example of such a 'defect' state 
is shown in Fig. [7| two localized magnons reside on a small \(fi) and a large \(Pl) 
hexagons. The large hexagon state is decomposed into a basis hexagon states (|2-5|) 
as 

6 

\<Pl) = + \<Pi) • (4-5) 
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Fig. 7. The defect state of two localized magnons, which breaks the hard-hexagon rule. 



The two-particle state {(fLVi) corresponds, therefore, to a special linear combination 
of states (|4-l(l and state \<pf) and, obviously, breaks the hard-hexagon rules. By- 
analogy with Fig. d a general defect two-magnon state \<pl<Pl') is constructed by 
drawing two closed lines L and L' on a kagome lattice such that (i) the lines can 
host localized magnons and (ii) one of the lines L' lies inside the area Al enclosed by 
the second line: L' C Al- Such states are, however, topologically equivalent to the 
two-magnon state in Fig. [7| and can be represented as linear combinations of normal 
hard-hexagon states and the primitive defect state. In addition, if one or both lines 
L and V have nonzero winding around the cluster, the corresponding two-magnon 
state cannot again be decomposed into hard-hexagon states (see discussion in Sec. 2). 
In the two-magnon subsector the total number of states breaking the hard-hexagon 
rule is oc N, that is a factor of 1/N smaller than the number of hard-hexagon states. 
Similar estimate holds for the three-magnon subsector: 0(N 2 ) defect states versus 
0(N 3 ) hard-hexagon states and so on. In addition, in magnetization subsectors with 
finite densities of magnons the number of defect states is further reduced. Indeed, the 
defect state shown in Fig.[7|exists only if there are no localized magnons on any of 12 
small hexagons surrounding the larger loop. As the total density n/N of (localized) 
magnons grows, there is less and less space for the composite defect states. Finally, if 
there are more than two magnons on every 19 hexagons or n/N > 2/(19 • 3) ~ 0.035, 
the number of defect states becomes significantly suppressed and eventually goes to 
zero. 

Another complication for the kagome antiferromagnet comes from the fact that 
the gap for one of the dispersive branches in Eq. (|2-2() also vanishes at the saturation 
field. The corresponding propagating magnon has the same energy ^(O) = H—3J as 
localized magnons from the flat branch. The lowest-energy states in the two-magnon 
subsector contain apart from the isolated localized magnons also superposition states 
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of one localized magnon and one propagating magnon. Such states form a continuum 
above the low-energy threshold Ei = 2 (if — 3J). The same is true for all n-magnon 
subsectors with n <^ N. Once n becomes a finite fraction of iV the above picture 
changes. The low-energy propagating magnons experience multiple scattering from 
an infinite number of localized magnons. Such scattering produces a finite shift of the 
energy of propagating magnons. An exact value of energy shift for 0*3 (k) depends on 
a precise pattern of localized magnons. An explicit estimate can be obtained for an 
average 'uniform' state of n-localized magnons, corresponding to a low-temperature 
ensemble of states with different translational patterns. We bosonize the spin Hamil- 
tonian using, e.g., the Dyson-Maleev transformation: = y/2S(l — b\bi/2S)bi, 
= y/2Sb\, for an arbitrary value of on-site spin S. Diagonalization of quadratic 
terms gives the excitation spectrum Q2 2|) . Effect of interaction between magnons 
is treated in the lowest order by using the Hartree-Fock approximation. We define 
two averages: for an on-site density m\ = (b\bi) and for a nearest-neighbor hop- 
ping m-2 = (b\bj). After a mean-field decoupling of four-boson terms the effect of 
interaction is reduced to a renormalization S — ► (S — mi + 1712) in the quadratic 
terms. The lowest point of the dispersive branch is at 0^3(0) = H — 6J(S — mi + 1712), 
which reduces at the saturation field H s = 6JS to Ad = 6 J (mi — m.2). For a single 
localized magnon (|2-5[) one finds mi = — 7722 = 1/6. Hence, the dispersive mode 
propagating through an 'averaged ensemble' of n localized magnons acquires a fi- 
nite gap Ad = 3Jn/N, which separates the dispersive branch from the localized 
magnons. At low enough temperatures T <C Ad the dispersive modes have a negli- 
gible contribution to the thermodynamics of a kagome antiferromagnet. Note, that 
characteristic temperature, where the above approximation becomes valid, depends 
on an average density of magnons and, therefore, on an applied magnetic field. 

4.2. Low-temperature behavior 

In the above analysis we have established that the hard-hexagon description of 
lowest energy states is valid at sufficiently high densities of magnons. High or low 
means here in comparison to the density of the magnon crystal: n/N = 1/9. Let us 
now discuss the thermodynamic properties of a spin-1/2 kagome antiferromagnet in 
the hard-hexagon approximation. The partition function of the hard-hexagon lattice 
gas is given by the same expression (|3-5|) . where the site index runs over J\f = N/3 
sites of a triangular lattice formed by centers of hexagon voids of the original kagome 
lattice. This model was solved by Baxter with the help of a corner transfer matrix 
method. 18 ^ Specifically, the normalized partition function Z 1 ^ and the fugacity 
z = e^/ T are expressed as functions of a real auxiliary parameter x: at the high-field 
region < z < z c (— 1 < x < 0), 



Z = - X %M> H(x) = H[(l-x 5n - i )(l-x 5n - l )]-\ (4-6) 

^ ' 71=1 

OO OO 

G(x) = I] [(1 - x 5 "" 3 )(l - X 5 "" 2 )]- 1 , Q(x) = H(l-x n ), 

71=1 71=1 
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l/N _ H 3 {x)Q 2 {x 5 ) fr (1 - x 6 "~ 4 )(l - x 6n ~ 3 ) 2 (l - x 6n ~ 2 ) 
G 2 (x) 11 (l-x 6 ™- 5 )(l-x6«-i)(l-x 6 ™) 2 ; 

v ' n=l ' 

whereas at the low-field region z > z c (0 < x < 1), 

flfx) - Q( " )Q(X5) (4-7) 
Z " X H*(x) ' 1 ' ~ Q 2 (x3) ' 14 '> 



m x-^G\x)Q\x^) fr (l-x 3 "- 2 )(l 

H 2 (x) J-- 1 - (1-x 3 ™) 2 

n=l ' 



A nonanalytic behavior of the partition function at z = z c = (ll+5\/5)/2 « 11.09017 
as x — ► ±1 indicates presence of a critical point. Such a continuous order-disorder 
transition corresponds to a spontaneous occupation of one of the three triangular 
sublattices at high densities. Expansion of Z in powers of z — z c yields exact critical 
exponents of the hard-hexagon model, 18 ) whereas numerical evaluation of rapidly 
converging infinite products provides dependence of various physical quantities on 
the chemical potential (temperature). 

Low-temperature behavior of a spin- 1/2 kagome in the vicinity of the saturation 
field can be obtained from the above exact solution by a trivial rescaling from the 
number of hexagons to the number of spins and using the relation \x = H s — H. The 
entropy as a function of magnetic field at constant T reaches a sharp maximum at 
H = H s . At the saturation field (z = 1, x = —0.25496) both the entropy and the 
magnetization have universal temperature independent values: 

S s /N = 0.11108 w 0.16026 In 2 , M s /N = 0.44596 . (4-8) 

The above value of the magnetization corresponds to the density n s /N ~ 0.054 
of localized magnons. This density well exceeds an estimate n ~ 0.035 from the 
previous subsection, which signifies suppression of the non-hard-hexagon states. As 
an external field is further decreased there is a continuous crystallization transition 
into a close-packed structure (Fig. |5J) with broken translational symmetry. The 
transition takes place at a temperature dependent critical field H C (T): 

H C (T) = H s -Tlnz c « 3 J - 2.40606T , (4-9) 

whereas the magnetization has again a universal value: 

M c /N = 1/2- 1/3 p c , Pc = (5 - VE)/10 . (4-10) 

The density of magnons at the transition point n c /N = p c /3 ~ 0.0921 is about 20% 
smaller than the density of the ideal close-packed structure n/N =1/9. The second 
order transition at H C (T) belongs to the universality class of a two-dimensional 
three-state Potts model 19 ) and has the following exact critical exponents: a = 1/3, 
(3 = 1/9, 7 = 13/9, and v = 5/6. 18 ^ The uniform susceptibility x = dM/dH 
also diverges at the transition point with the same critical exponent a as the specific 
heat. Apart from a sharp lambda-peak at z = z c , the specific heat has also a rounded 
Schottky type anomaly at z m = 0.03897 (x = -0.03296, lnz m = -3.24502). 
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Fig. 8. Field dependence of the specific heat (left) and the magnetization (right) of a spin-1/2 
kagome antiferromagnet obtained in the hard-hexagon model. 

Hard-hexagon results for the field dependence of the specific heat and the mag- 
netization of a spin-1/2 kagome antiferromagnet are shown in Fig. |HJ Positions of 
two peaks in the specific heat scale linearly with temperature. As a function of 
temperature at a fixed magnetic field the specific heat also exhibits a characteristic 
low-temperature peak at T m = (H s — H)/lnz m , see Fig. [SJ At higher temperatures 
T ~ J there should be a second maximum in C(T), which is not captured by hard- 
hexagon mapping and corresponds to development of short-range spin correlations. 
Overall, the behavior of the specific heat and the magnetization of a kagome antifer- 
romagnet resembles the previous results for the sawtooth chain. The major difference 
is that while one-dimensional model has only a finite temperature crossover, in two- 
dimensions the breaking of a discrete translational symmetry occurs via a sharp 
second-order transition. Experimental observations of a double peak structure in 
a temperature (field) dependence of the specific heat and of a characteristic non- 
monotonic temperature dependence of the magnetization (Fig. |UJ right panel) could 
be a clear signature of degenerate magnons in a geometrically frustrated magnetic 
material. 

Since the gaps between the noninteracting localized magnons and propagating 
states are smaller for a kagome antiferromagnet, the quantitative validity of the hard- 
hexagon mapping is restricted to lower temperatures T < 0.05-0.1J and to higher 
magnetic fields H > 0.9ff s . The critical behavior obtained from the exact solution 
for hard hexagons should, however, remain valid even beyond the above range. At 
present there are no numerical results on the magnetothermodynamics of a spin-1/2 
kagome antiferromagnet in order to compare them with the hard-hexagon mapping. 
Note, that exact diagonalization results on finite periodic clusters should suffer from 
rather strong finite size effects. They are determined by a presence of large number 
of non-hard-hexagon states of localized magnons, such as the 'defect' states discussed 
earlier in this section, whose contribution vanishes only in the thermodynamic limit 
N -> oo. 

The order parameter of the hard-hexagon model is defined as a difference of 
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Fig. 9. Temperature dependence of the specific heat (left) and the magnetization (right) of a spin- 
1/2 kagome antiferromagnet obtained in the hard-hexagon model. Curves on the left panel 
correspond to H — 3.1J, 3.25J, 3.5J, and 4J from left to right. 

hexagon densities pi on adjacent sites: R = pi — pi +a . 18 ^ Up to a renormalization 
prefactor R is equal to the Fourier harmonics p q at q = (47r/3a*,0), where a* is a 
period of a triangular lattice. In a kagome antiferromagnet the magnon crystal has 
two types of spins, on hexagons with localized magnons and between them. The 
two sublattices have different average magnetizations: (S z ) = 1/3 and 1/2 at T = 0. 
Hence, the crystalline order of magnons reveals itself as an extra Fourier harmonics 
in the magnetic structure factor S zz (q). The wave- vector of magnon crystal is the 
same one as for a so-called \/3 x \/3 structure, stabilized in zero field for a kagome 
antiferromagnet with large S 3> l. 20 ) Note, that the magnon crystal state is not 
reduced to a semiclassical collinear two-sublattice order. The quantum coherence 
of spin-flip propagating around one hexagon plays a crucial role in stabilizing the 
magnon crystal. 

Finally, the schematic phase diagram of a Heisenberg spin-1/2 kagome antifer- 
romagnet in external magnetic field is shown in Fig. I1UI In zero magnetic field 
this quantum frustrated model has a nonmagnetic ground state with a large num- 
ber of low-lying singlet excitations between the ground state and the lowest excited 
triplet. 21 )' 22 ) Presence or absence of a finite-temperature phase transition at H = 
depends on whether such a singlet ground state is a valence bond crystal, which 
breaks certain discrete lattice symmetries, or a true spin-liquid state with full sym- 
metry of the Heisenberg Hamiltonian. The answer to this question is not known 
at present due to a limited size of clusters accessible for numerical diagonalization. 
Applied magnetic field eventually closes the triplet gap and the system transforms 
into a state with finite magnetization. Numerical investigations of the magnetization 
curve of a spin-1/2 Heisenberg model find a clear plateau at 1/3 of the saturation 
magnetization. 23 )' 24 - 1 The ground state at the 1/3-plateau has a long-range order 
of valence-bond type with the same \/3 x \/3 periodicity as the magnon crystal in 
Fig. El 25 -* There is, therefore, a line of phase transitions, which separates the 1/3 
plateau state from a paramagnetic phase and phases above and below the plateau. 
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Fig. 10. Schematic H-T phase diagram of a Heisenberg spin-1/2 kagome antiferromagnet. Dashed 
area shows a region of applicability of the hard-hexagon mapping. 

The magnon crystal state is the third presently known phase of a spin-1/2 kagome 
antiferromagnet, which is stabilized in the vicinity of the saturation field. At zero 
temperature the magnon crystal corresponds to the 7/9-magnetization plateau. The 
hard-hexagon mapping obtained in Ref. 7) and further developed in the present work 
is valid in the dashed region of the phase diagram. Being restricted to high fields 
and low temperatures the hard-hexagon mapping is, nevertheless, quite important, 
since it provides the only exact information about the phase diagram of a spin-1/2 
Heisenberg antiferromagnet on a kagome lattice. 

Between the above three phases, magnetization of a kagome antiferromagnet 
grows continuously, 23 -*' 24 ) which points at a gapless excitation spectrum. The number 
and the nature of such intermediate phases remain an open issue. Gapless states with 
finite magnetization can also have nontrivial properties and exotic order parameters 
even in the classical limit S — > oo. 26 ) 

§5. Discussion 

The above analysis of the magnetothermodynamics of the sawtooth chain and 
of the quantum kagome antiferromagnet has been focused mainly on spin-1/2 mod- 
els. Apart from a detailed presentation of the previous results 4 )' 6 -*' 7 ) we have have 
developed an effective Ising description of localized magnons for the sawtoth chain, 
which is valid up to significantly higher temperatures than the hard-dimer approx- 
imation. Most of the discussed results remain valid for an arbitrary quantum spin 
5*. These include the universal values of the entropy at H = H s and the finite tem- 
perature transition into the magnon crystal state for a kagome lattice model. The 
zero temperature magnetization jumps become, however, increasingly small for large 
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S: a relative height of the jump decreases as AM = 1/(95) for a kagome antifer- 
romagnet. The temperature range for stability of the magnon crystal also becomes 
smaller with increasing S. Indeed, the gap between localized magnon states and 
higher-energy propagating states is proportional to the exchange constant J. (Vari- 
ational calculations for the S = 1 sawtooth chain yield only a slightly different gap 
A = 0.41J compared to the S = 1/2 result of Sec. 3.) The crystallization tempera- 
ture scales, therefore, as T c = O(J) and becomes a small fraction of the Curie- Weiss 
constant T c /9cw ~ 1/S 2 . The domain of quantum effects is pushed to extremely 
low temperatures for semiclassical spins S S> 1 

A strong magneto caloric effect discussed in Sec. 3 remains, however, present 
even in the classical limit S — ► oo. 16 ) Recently, a large temperature decrease during 
an adiabatic demagnetization process has been measured in pyrochlore compound 
Gd2Ti2C>7 with S = 7/2. 27 ) Note, that antiferromagnets on a pyrochlore and a 
checkerboard lattices also have localized magnon excitations. Their low-temperature 
behavior is, therefore, mapped to lattice gas models of hard-core classical particles. 
For the checkerboard antiferromagnet an appropriate low-temperature model is a 
square lattice gas of hard-core particles with nearest- and next-nearest neighbor 
exclusions. There is no exact solution for such statistical model. Still, effective 
representations may be quite useful, because lattice gas models can be efficiently 
studied with a powerful classical Monte Carlo method. The thermodynamics of 
the corresponding lattice models generally depends on a dimensionless fugacity z = 
e (H 3 —H)/T _ ^his immediately implies a strong magneto caloric effect T — > as H —* 
H s for all frustrated magnets with localized spin-flips. The residual entropy at 
H = H s , which determines a cooling power, does depend on a lattice geometry. 
An interesting open problem is calculation of the residual entropies for pyrochlore 
and checkerboard antiferromagnets in order to determine the favorable geometrically 
frustrated spin system for applications in adiabatic refrigerators. 28 ^' 29 ^ There is also 
a certain difference in the magnetocaloric effect for quantum and (semi) classical 
frustrated magnets. The quantum spin systems exhibit a strong cooling on the two 
sides of the saturation field when H — > H s db 0: both the saturated phase at H > H s 
and the magnon crystal at H < H s have no entropy at T = 0. Frustrated classical 
spin system has in contrast infinite degeneracy of the ground state for H < H s . 
Therefore, a much stronger cooling effect is found when applied field is decreased 
towards H s from the high-field side. 16 ) For intermediate values of spin S > 1 the 
quantum order by disorder effect 30 ) lifts the classical degeneracy at H < H s and 
restores zero entropy. This leads to an asymmetric field dependence for Tg(H) with 
the steepest slope at H > H s and the lowest temperature achievable in the vicinity 
of the saturation field. 

Real magnetic materials usually have additional weak interactions apart from a 
nearest-neighbor Heisenberg exchange. These include, for example, various anisot- 
ropics, next-nearest-neighbor exchanges, and the dipolar interactions. Such extra 
interactions lift the classical degeneracy of frustrated magnets and in the quan- 
tum case induce a weak dispersion of the lowest magnon branches ()2-2|) and (|2-4|) . 
Phase transition at H = H s corresponds, then, to a usual condensation of a single 
magnon mode at a certain wave- vector Q. The finite zero-temperature entropy at 
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the saturation field is completely destroyed in such a case. The question, whether 
some of the effects described in the present paper are still observable or not, can be 
solved by a simple comparison of the energy scales. If the gap in a magnon crystal 
state exceeds the bandwidth W of the lowest branch cJi(k), then there has to be a 
subsequent phase transition between a state with transverse magnetic order in the 
immediate vicinity of H s and the crystal state without such ordering. (There will 
be, of course, zero-point fluctuations in the magnon crystal ground state on top of a 
simple product of localized magnon wave- functions, Fig. |2J) As long as the energy 
scale for weak residual interactions is significantly smaller than the nearest-neighbor 
exchange W <C J, a frustrated magnet continues to exhibit a large magneto caloric 
effect. The lowest reachable temperature is, however, limited to 

-^min rsJ 

W. Finally, 

one of the suitable ways to detect experimentally nearly dispersionless magnon bands 
in frustrated magnets is to measure a low-temperature peak in the specific heat for 
T m = (H-H s )/lnz m >W. 
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